tableOlogit <- function(formulas,y=data,modelnames,...){
  require(MASS)
  foo <- vector("list",length(formulas))
  for(i in 1:length(formulas)) {
    foo[[i]] <- polr(substitute(i,list(i=formulas[i])), 
                     data=y,Hess=TRUE)
   }
  return(foo)
}
